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PsJ ' Abstract. A new method, based on simulated annealing technique and 

pH . aimed at the inverse problem in the analysis of intergalactic or interstellar 

I complex spectra of hydrogen and metal lines, is outlined. We consider 

• the process of line formation in clumpy stochastic media accounting for 

'nI" I fluctuating velocity and density fields self-consistently. Two examples 

■ of the analysis of 'H+D'-like absorptions seen at Zg, = 3.514 and 3.378 
^ ; towards APM 08279+5255 are presented. 

> ■ 

^ : 

O . 

. 1. Introduction 

o : 

. The cosmological significance of the deuterium abundance measurements in 

^ I metal-deficient QSO absorption systems has been widely discussed in the lit- 

■ erature (see e.g. the review by Lemoine et al. 1999). Practical applications 
Qh| of such measurements were clearly outlined in Tytler & Buries (1997) : (i) the 
Q ■ primordial D/H value gives the density of baryons Sib at the time of big bang 
|H I nucleosynthesis (BBN), a precise value for might be used (ii) to determine 
t/^ ■ the fraction of baryons which are missing, (iii) to specify Galactic chemical evo- 

\ lution, and (iv) to test models of high energy physics. A measurement of the 

^ ■ D/H ratio, together with other three light element abundances (^He/H, ^He/H, 

^ I and '^Li/H), provides the complete test of the standard BBN model. 

^ • Deuterium has been reported up-to-now in a few QSO Lyman limit systems 

, (LLS), - the systems with the neutral hydrogen column densities of A'^hi ~ 10^^ — 

10^*^ cm~^ (Kirkman et al. 1999). The difficulties inherent to measurements of 
D/H in QSO spectra are mainly caused by the confusion between the Di line 
(which is always partially blended by the blue wing of the saturated hydrogen 
line) and the numerous neighboring weak lines of H i observed in the Lya forest 
at redshifts z>2 (Buries & Tytler 1998). 

Currently, there are two methods to analyse the absorption spectra : (i) a 
conventional Voigt-profile fitting (VPF) procedure which usually assumes sev- 
eral subcomponents with their own physical parameters to describe a complex 
absorption profile, and (ii) a mesoturbulent approach which describes the line 
formation process in a continuous medium with fluctuating physical characteris- 



^ Based on data obtained at the W. M. Keck Observatory, which is jointly operated by the 
University of Cahfornia, the California Institute of Technology, and the National Aeronautics 
and Space Administration. 



1 



2 



Levshakov S. A. 



tics. It is hard to favor this or that method if both of them provide good fitting. 
But the observed increasing complexity of the Une profiles with increasing spec- 
tral resolution gives some preference to the model of the fluctuating continuous 
medium. 

Here, we set forward a mesoturbulent approach to measure D/H and metal 
abundances, which has many advantages over the standard VPF procedures. 
A brief description of a new Monte Carlo inversion (MCI) method is given in 
this report. For more details, the reader is referred to Levshakov, Agafonova, &; 
Kegel (2000b). 

An example of the MCI analysis of two 'H+D'-like profiles with accompa- 
nying metal lines observed at Za = 3.514 and 3.378 towards the quasar APM 
08279+5255 is described. The high quality spectral data have been obtained 
with the Keck-I telescope and the HIRES spectrograph by Ellison et al. (1999). 

2. The MCI method and results 

The MCI method is based on simulated annealing technique and aimed at the 
evaluation both the physical parameters of the gas cloud and the corresponding 
velocity and density distributions along the line of sight. We consider the line 
formation process in clumpy stochastic media with fluctuating velocity and den- 
sity fields (mesoturbulence) . The new approach generalizes our previous Reverse 
Monte Carlo (Levshakov, Kegel, & Takahara 1999) and Entropy-Regularized 
Minimization (Levshakov, Takahara, &; Agafonova 1999) methods dealing with 
incompressible turbulence (i.e. the case of random bulk motions with homoge- 
neous gas density nn and kinetic temperature T). 

The main goal is to solve the inverse problem, i.e. the problem to deduce 
physical parameters from a QSO absorption system. The inversion is always an 
optimization problem in which an objective function is minimized. To estimate 
a goodness of the minimization we used a function augmented by a regular- 
ization term (a penalty function) to stabilize the MCI solutions. The MCI is 
a stochastic optimization procedure and one docs not know in advance if the 
global minimum of the objective function is reached in a single run. Therefore 
to check the convergency, several runs are executed for a given data set with 
every calculation starting from a random point in the simulation parameter box 
and from completely random configurations of the velocity and density fields. 
After these runs, the median estimation of the model parameters is performed. 

Our model supposes a continuous absorbing gas slab of a thickness L. The 
velocity component along a given line of sight is described by a random function 
in which the velocities in neighboring volume elements are correlated with each 
other. The gas is optically thin in the Lyman continuum. We are considering 
a compressible gas, i.e. nn is also a random function of the space coordinate, 
X. Following Donahue & Shull (1991) and assuming that the ionizing radiation 
field is constant, the ionization of different elements can be described by one 
parameter only - the ionization parameter U oc l/nn- Furthermore, for gas in 
thermal equilibrium, Donahue &; Shull give an explicit relation between U and 
T. The background ionizing spectrum is taken from Mathews Sz Ferland (1987). 

In our computations, the continuous random functions v{x) and the nor- 
malized density y{x) = nii{x)/no , uq being the mean hydrogen density, are 



D/H in QSOs : a mesoturbulent approach 



3 



represented by their sampled values at equally spaced intervals Ax, i.e. by the 
vectors {vi, ■.■,Vk} and {yi, yfc} with k large enough to describe the narrowest 
components of complex spectral lines. For the ionization parameter as a func- 
tion of X, we have U{x) = Uo/y{x), with Uq being the reduced mean ionization 
parameter defined below. 

Absorption system at Za = 3.514. A measurement of the primordial D/H in 
the Za, = 3.514 system has been recently made by Molaro et al. (1999). They 
suggested that the blue wing of H l Lya is contaminated by D I and evaluated 
a very low deuterium abundance of D/H~ 1.5 x 10"^ in the cloud with A'^hi = 
(1.23lQ'Qg) X 10^^ cm~^. They considered, however, the derived D abundance 
as a lower limit because their analysis was based on a simplified one-component 
VPF model which failed to fit the red wing of the Lya line as well as the profiles 
of Sim, Si IV, and Civ lines exhibiting complex structures over approximately 
100 km velocity range. They further assumed that additional components 
would decrease the H l column density for the major component and, thus, would 
yield a higher deuterium abundance. Given the MCI method, we can test this 
assumption since the MCI accounts self-consistently for the velocity and density 
fluctuations. 

Our aim is to fit the model spectra simultaneously to the observed Hi, C II, 
Civ, Sim, and SilV profiles. In this case the mesoturbulent model requires the 
definition of a simulation box for the six parameters : the carbon and silicon 
abundances, Zq and Z^i, respectively, the rms velocity and density dispersion 
CTy, the reduced total hydrogen column density Ah = Ah/(1 + Cy)^/^, and the 
reduced mean ionization parameter Uq = C/o/(l + o-y)^''^- For the model param- 
eters the following boundaries were adopted : Zc ranges from 10-^ to 4 X 10"^, 
Zsi from 10"^ to 3 X 10-^ ay from 25 to 80 km s'^, ay from 0.5 to 2.2, iVn 
from 5 X 10^^ to 8 x 10^^ cm , and Uo ranges from 5 X 10"^ to 5 x lO'^. We 
fix Zg, = 3.51374 (the value adopted by Molaro et al.) as a more or less arbitrary 
reference velocity at which Vj = 0. 

Having specified the parameter space, we minimize the value. The ob- 
jective function includes those pixels which are critical to the fit. In Fig. 1, 
these pixels are marked by shaded areas. In Fig. 1 (panels f and c), the ob- 
served profiles of C ivA1550 and, respectively, Si ivA1402 are shown together with 
the model spectra computed with the parameters derived from Lya, C iiA1334, 
CivA1548, SiiiiA1206, and SiivA1393 fitting to illustrate the consistency. For 
the same reason the Ly/3 model spectrum is shown in panel b at the expected 
position. All model spectra in Fig. 1 are drawn by continuous curves, whereas 
filled circles represent observations (normalized fluxes). The corresponding dis- 
tributions of v{x), y{x), and T{x) are shown in panels i, j, and k. The restored 
velocity field reveals a complex structure which is manifested in non-Gaussian 
density-weighted velocity distribution as shown in panels 1 and m for the total 
hydrogen as well as for the individual ions. We found that the radial velocity 
distribution of Hi in the vicinity of Av ~ —100 km s~^ may mimic the deu- 
terium absorption and, thus, the asymmetric blue wing of the hydrogen Lya 
absorption may be readily explained by H l alone. 

The median estimation of the model parameters gives A^h = 5.9 x 10^^ 
cm-2, Ahi = 5.3 X 10^5 cm'^, Uq = 1.6 x lO'^, ay = 51 km s'^, and ay = 1.1. 
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Figure 1. (a-h). Observed normalized intensities (dots) and theo- 
retical (solid curves) profiles of the absorption lines from the z^, = 3.514 
system (note different scales in panels h and f). Shaded areas mark 
pixels which are critical to the MCI fit. (i-k). The corresponding MCI 
reconstraction of the radial velocity, density, and kinetic temperature 
T4 which is given in units of 10^ K. (l,m). The corresponding density- 
weighted radial velocity distribution functions, p{v), for the total Htot 
and neutral Hi hydrogen, Cll, CiV, Sim, and Siiv as restored by the 
MCI. For comparison, in panel 1 the short dashed curve shows p{v) 
adopted by Molaro et al. 
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Figure 2. (a-d). Observed normalized intensities (dots) and theo- 
retical (solid curves) profiles of the absorption lines from the Za = 3.378 
system (shaded areas mark pixels used by the MCI). The 'H-fD'-like 
absorption is explained by H I alone. An excellent common fit to H I, 
Si IV, and Civ should, however, be rejected since the upper panel 
demonstrates the inconsistency between the synthetic spectrum of C ii 
computed with the parameters derived from the foregoing lines and the 
observed intensities. 

The results were obtained with k = 100 and the correlation coefficients /y = 
/y = 0.95 (for more details, see Levshakov, Agafonova, & Kegel 2000a). 

The MCI allowed us to fit precisely not only the blue wing of the saturated 
Lya line but the red one as well. We found that the actual neutral hydrogen 
column density may be a factor of 250 lower than the value obtained by Molaro 
et al. if one accounts for the velocity field structure. Besides we did not confirm 
the extremely low mctallicity of [C/H] ~ —4.0, and [Si/H] ~ —3.5 reported by 
Molaro et al. Our analysis yields [C/H] ~ —1.8, and [Si/H] ~ —0.7. A similar 
silicon overabundance has also been observed in halo (population II) stars (Henry 
& Worthey 1999). 

Absorption system at Zg, = 3.378. The following example illustrates how the 
realibility of the inversion procedure can be controlled. We have chosen the 
Za = 3.378 system since at the position of the narrowest CiV subcomponent 
with zqiv = 3.37757 and 6c iv = 6.5 km s^^ (see Ellison et al.) one can see an 
'H+D'-likc absorption in the blue wing of the saturated Lya line (Fig. 2a). The 
C IV, and Si IV profiles from this system were treated by Ellison et al. separately. 
They found Ar^°V = i.^-'^^-oii) x 10^^ cm^^ and iV^?V = (1-70 ± 0.08) x 10^^ 
cm~^. For the neutral hydrogen, they estimated A^hi = 2.8 x 10^^ cm~^ and 
6hi = 78.6 km s~^ (errors for both quantities are greater than 30%). 

In this example, we assumed no deuterium absorption and tried to force a 
common fit to the lines shown in panels a-c (Fig. 2). Pixels used in the fitting 
are labeled by shaded areas. The MCI fit, shown by solid curves, looks excellent, 
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and gives TVhi = 2.1 x 10^^ cm-^, iVdv = 1.0 x 10^^ cm-^, A^gij^ = 2.6 x lO" 
cm~^, cjv = 41.6 km s^^, and Cy = 1.3. The results were obtained with k = 150 
and the correlation coefficients /v = /y = 0.97. 

The obtained MCI solution should, however, be rejected because the syn- 
thetic spectrum of Cii (solid curve in Fig. 2d), computed with the parameters 
derived from the Lya, Civ, and Siiv fitting, differs significantly from the ob- 
served intensities (dots in Fig. 2d). This example shows that we can always 
control the MCI results using additional portions of the analysed spectrum. 

Another issue of this example is that we can, in principle, fit an 'H+D'-like 
absorption by Hi alone even for the systems with A^hi — 10^^ cm~^ and accom- 
panying metal lines. Examples of false deuterium identifications in systems with 
10 times lower neutral hydrogen column densities and without supporting metal 
lines have been discussed in Tytler & Buries. Both cases stress the importance 
of the comprehensive approach to the analysis of each individual QSO system 
showing possible D absorption. 

We may conclude that up-to-now deuterium was detected in only four QSO 
spectra (Q 1937-1009, Q 1009+2956, Q 0130-4021, and Q 1718+4807) where 
Ah I was measured with a sufficiently high accuracy. These measurements are 
in concordance with D/H = (3 — 4) x 10~^. 
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